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Abstract 

The Yang-Lee zeros of the Q-state Potts model on recursive lattices are studied for non- 
integer values of Q. Considering ID lattice as a Bethe lattice with coordination number 
equal to two, the location of Yang-Lee zeros of ID ferromagnetic and antiferromagnetic 
Potts models is completely analyzed in terms of neutral periodical points. Three different 
regimes for Yang-Lee zeros are found for Q > 1 and < Q < 1. An exact analytical 
formula for the equation of phase transition points is derived for the ID case. It is shown 
that Yang-Lee zeros of the Q-state Potts model on a Bethe lattice are located on arcs of 
circles with the radius depending on Q and temperature for Q > 1. Complex magnetic 
field metastability regions are studied for the Q > 1 and < Q < 1 cases. The Yang-Lee 
edge singularity exponents are calculated for both ID and Bethe lattice Potts models. 
The dynamics of metastability regions for different values of Q is studied numerically. 
PACS number(s): 05.50.+q, 05.10.-a 

1 Introduction 

The Q-state Potts model plays an important role in the general theory of phase transitions 
and critical phenomena (![]. It was initially denned for an integer Q as a generalization of 
the Ising model (Q = 2) to more-than-two components [f2J. Later on, it was shown that the 
Potts model for non-integer values of Q may describe the properties of a number of physical 
systems such as dilute spin glasses 0, gelatation and vulcanization of branched polymers 
(0 < Q < 1) Also it was shown that the bond and site percolation problems could be 
formulated in terms of Potts models with pair and multisite interactions in the Q = 1 limit. 

In 1952, for the first time, Lee and Yang |j] in their famous papers studied the distribution 
of zeros of the partition function considered as a function of a complex magnetic filed (e~ 
activity, if is a magnetic field). They proved the circle theorem which states that zeros of 
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the partition function of an Ising ferromagnet lie on unit circle in the complex activity plane 
(Yang-Lee zeros). After these pioneer works of Lee and Yang, Fisher ||, in 1964, initiated 
the study of partition function zeros in the complex temperature plane (Fisher zeros). These 
methods were then extended to other types of interactions and found wide applications @. 

Recently, much attention was drawn to the problem of Yang-Lee and Fisher zeros of the 
Q-state Potts model for both integer and non-integer values of Q. The microcanonical transfer 
matrix method was used to study the Yang-Lee and Fisher zeros of the non-integer Q-state 



Potts model in two and three dimensions M- 11 



Derrida, De Seze and Itzykson showed for the first time that the Fisher zeros in hierar- 
chical lattice models are just the Julia set corresponding to the renormalization transformation. 
They found a fractal structure for the Fisher zeros in Q-state Potts model on the diamond lat- 
tice. Recently, Monroe investigated Julia sets of the Potts model on recursive lattices [p~3j| . 
He found that the box counting fractal dimension of the Julia set of the governing recurrence 
relation is a minimum at a phase transition. This gives an alternative criterion for studying 
phase transitions in models defined on recursive lattices. 

In 1994, Glumac and Uzelac WM, using the transfer matrix method, analytically studied 



the distribution of Yang-Lee zeros for the one-dimensional ferromagnetic Potts model with 
arbitrary and continuous Q > 0. For < Q < 1 they obtained that for high temperatures the 
Yang-Lee zeros lie on a real interval and for low temperatures these are located partially on the 
real axis and in complex conjugate pairs on the activity plane. Later on, Monroe investigated 
this model by means of the dynamical systems approach and confirmed that for Q < 1 there 
is a real interval of Yang-Lee zeros ||15|| . Then, Kim and Creswick found that for Q > 1 the 
Yang-Lee zeros lie on a circle with radius R, where R < 1 for l<Q<2,i?>lforQ>2 and 
R = 2 for Q — 2 |J. However, it is not clear yet what is the location of Yang-Lee zeros for 
< Q < 1 at low temperatures flUf] . 

In this paper the Yang-Lee zeros of the one-dimensional and Bethe lattice Q-state Potts mod- 
els are studied using the dynamical systems approach. It is shown that for the one-dimensional 
Potts model the partition function becomes zero when the corresponding recurrence relation 
has neutral fixed points for a given value of magnetic field. Using this correspondence between 
zeros of the partition function and neutral fixed points of the recurrence relation the Yang-Lee 
zeros of both ferromagnetic and antiferromagnetic Potts models are completely studied analyt- 
ically. The location of Yang-Lee zeros of the ferromagnetic Potts model for < Q < 1 is found. 
Also, formulas for the density of Yang-Lee zeros are derived and edge singularity exponents 
are calculated. For the Potts model on a Bethe lattice it is shown that the Yang-Lee zeros are 
located on a phase coexistence line in the complex magnetic field plane. Here, the phase coex- 
istence line is defined as a line in the complex magnetic field plane, where the absolute values of 
derivatives of the recurrence relation in two attracting fixed points are equal. It is noteworthy 
that, Monroe fl6j used a similar criterion for studying critical properties of the Potts model on 
recursive lattices. For the Bethe lattice case an analytical study of Yang-Lee zeros is performed 
also. A formula for Yang-Lee edge singularity points is derived and edge singularity exponents 
are calculated. Our analitical treatment confirmed the results of numerical calculations. Fur- 
ther, metastability regions in a complex magnetic field plane are investigated. It is shown that 
the border of a metastability region may be found from the condition of existence of a neutral 
fixed point of the governing recurrence relation. 

The structure of this paper is as follows, in Section 2 an exact recurrence relation (Potts- 
Bethe mapping) for the Q-state Potts model on a Bethe lattice is derived. Applying the theory 
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Figure 1: The Bethe lattice with coordination number 7 = 3. 

of dynamical systems to the problem of phase transitions it is shown that critical points may 
be associated with neutral periodical points of the corresponding mapping. In section 3 the 
Yang-Lee zeros and edge singularities of ferromagnetic and antiferromagnetic Potts models 
are studied analytically for non-integer Q, considering a ID lattice as a Bethe lattice with 
coordination number 7 = 2. In Section 4 the Yang-Lee zeros and edge singularities of the Potts 
model on a Bethe lattice with coordination number 7 > 2 are studied numerically In the last 
section the dynamics of complex magnetic metastability regions is studied numerically and the 
explanation of results is given. 

2 The Q-state Potts Model on the Bethe lattice 

The Q state Potts model in the magnetic field is defined by the Hamiltonian 

—(3H = J £ SfaaJ + h'ESfaO), (1) 

<i,j> i 

where oi takes the values 0, 1, 2, . . ., Q — 1 and (3 = 1/kT. The first sum in the r.h.s. of ([[]) 
goes over all edges and the second one over all sites on the lattice. The partition function of 
the model is given by 

{«} 

where the summation goes over all configurations of the system. 

Using the recursive structure of the Bethe lattice (Fig. |1]) one can derive an exact recur- 
rence relation and apply the theory of dynamical systems to investigate the thermodynamical 



3 



properties of models defined on it [T^J. Cutting the lattice at the central site one will ob- 
tain 7 interacting nth generation branches |18| . By denoting the partition function of the nth 
generation branch with the basic site in the state <To as g„(cr ) the partition function may be 
written as follows 

Z n =J2 exp{h5(a ,0)} \g n (a )}\ (2) 
W 

where <7o is the Potts variable at the central site of the lattice (Fig. [I]). Applying the "cutting" 
procedure to a nth generation branch one can derive the recurrence relation for g n (cro) 

9n(<To) = exp{J5(a ,ai) + h5(a u 0)} [g n - 1 {(T 1 )\ r '' . (3) 



Introducing the notation 
one can obtain the Potts-Bethe map from the Eq(||D 



x„ = ^|, (4) 
9n{cr = 0) 



ft \ ft \ /^+(^ + Q-2) x^ 1 

where \x = e h , and z = e J . 

The magnetization of the central site for a Bethe lattice of n generations may be written as 

^ = ^g%..0)«-" % + (o ''_ 1)ri . (6) 

Instead of M(//) in Eq.@ it is more convenient to consider the function M(/i) = 2M(fi) — 1, 
which has the same analytical properties as M(/i) and give correct magnetization for the Ising 
model (Q = 2) 

MM = ^g-f . (7) 

In the following the formulas in Eqs.@ and @ will be generalized to non-integer values of 
Q>0. 

Let us now give some definitions and briefly discuss the problem of phase transitions on 
recursive lattices in terms of dynamical systems theory. The point x* is called a periodical 
point with period n of the mapping x n = /(x n _i) if it is a solution to the equation f k (x) = x. 
Here f n means a superposition f n = fofo---of. If k = 1, x* is called a fixed point. To 
analyze the stability of a periodic point x* of period k, i.e. either iterations of f(x) tend to the 
periodical point or move away from it, one should compute the derivative of f k , A = (f h )'(x) 
(' = iz)' a ^ this point. A periodic point x* is: (1) attracting (stable) if |A| < 1; (2) repelling 
(unstable) if |A| > 1 and (3) neutral (indifferent) if |A| = 1. 

The thermodynamic properties of models defined on recursive lattices may be investigated 
by studying the dynamics of the corresponding recursive function. For the ferromagnetic Ising 
model, for example, this may be done as follows. For high temperatures (T > T c ) the recursive 
function Eq.fl^) has only one attracting fixed point corresponding to a stable paramagnetic 
state. For low temperatures (T < T c ) the recursive function Eq.(0) has two attracting fixed 
points. In the absence of a magnetic field these two attracting fixed points correspond to two 
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possible ferromagnetic states with opposite magnetizations. It is well known that for h = 
and T < T c the system undergoes a first order phase transition and the condition |Ai| = | A2 1 
is satisfied, where (Ai^l are derivatives of f(x) in these attracting fixed points. In the presence 
of a magnetic field the stable state corresponds to the fixed point with maximum |A| and the 
magnetization in this state has the same direction as the external magnetic field. The other 
fixed point corresponds to a metastable state, which may be achieved by a sudden reversal 
in the sign of the magnetic field. The boundary of the metastability region may be found 
from the condition that one of the fixed points becomes neutral. For more details about 



the dynamics of metastable states see the book by Chaikin and Lubensky |p6fl . The critical 
temperature corresponds to the values of the magnetic field and temperature when the fixed 
points of recursive function f(x) are neutral and repelling. The values of external parameters 
(temperature, magnetic field, etc.) at which the recursive function f(x) has a neutral periodical 
point of period k may be obtained from the following system of equations 




Thus, for the models defined on recursive lattices, critical points and the boundary of metasta- 
bility regions may be associated with neutral periodical points of the model's mapping. In the 
following section the Yang-Lee zeros of the ID Q-state Potts model are investigated in terms 
of neutral fixed points. 



3 The Yang-Lee zeros of ID Q-state Potts model 

A one-dimensional lattice may be considered as a particular case of the Bethe lattice with 
coordination number 7 = 2. In this case the Bethe- Potts mapping (||) becomes a Mobius 
transformation, i.e. a rational map of the form 

ax + b 

R{x) = ad — be ^ 0, 

cx + a 

where 

R(oc) = a/c, R(—d/c) = 00, 



if c 7^ 0, while R(oo) = 00 when c = 0. The dynamics of such maps is rather simple ||19|| : If 
f(x) is a Mobius transformation with two fixed points, then either f n (x) converge to one of 
the fixed points of f(x) (one of the fixed points is attracting and the other is repelling), or 
they move cyclically through a finite set of points, or they form a dense subset of some circle 
(both fixed points are neutral). It follows from the dynamics of the Mobius transformation 
that phase transitions correspond to values of the temperature and magnetic field where the 
recursive function f(x) has only neutral fixed points. Thus, to study the Yang-Lee zeros one 
should study the system (H) to find the values of fi at which this system has solutions for k — 1. 
For neutral fixed points the system (^]) may be written in the form 

= e** 0g[O,2tt], ^ 
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Excluding x from the equations of the system (||) after some algebra one can find the equation 
of phase transitions in the form 

- 2 \(z - 1) + Q - 1) cos0 + 1 - Q] fj, + + Q - 2) 2 = 0. (10) 

where G [0, 27r]. Since cos0 is an even function of 0, we may restrict ourselves to G [0, 7r]. 
Thus, for given 2 and Q the equation fllCf ) is a parametric equation of Yang-Lee zeros, where 
is a parameter. Analyzing Eq. (|T0| ) one can find the location of Yang-Lee zeros of the Potts 
model. It is noteworthy that using this equation one can study the Fisher and Potts zerosQ in 
the same way[|. 

First of all note that Eq. fllOD is a quadratic equation of /x with real coefficients. Hence, 
solutions to this equation lie either on the real axes or in complex conjugate pairs on a circle 
with radius R = |/x| = l z+ Q~ 2 l for any G [0, ir] (See also 0). The solutions to the Eq (fL0|) can 
be written in the form 



= A [2 cos 2 § - B ± 2 ^cos 2 § (cos 2 f - 5)] , (11 



where 

A= («-l)(* + Q-l) and fl ^ + Q- 2 ) - (12) 

* 2 («-i)(z + g-i) 1 ; 

Let us now study the Yang-Lee zeros of ferromagnetic Potts model, i.e. J > and z > 1. 
For Q > 1 one can easily find that B > 1, hence, all solutions (O) are complex conjugate and 
lie on an arc of circle with radius i? = z+( ^~ 2 . If 1 < Q < 2 then i? < 1 (the arc lies inside 
the unit circle), if Q > 2 then _R > 1 (the arc lies outside the unit circle 0) and R = 1 |5| for 
Q = 2 (Ising model). Writing /x in the exponential form fi = Re 10 one can find 





cos ■ 



l)(z + Q-l) . 

cos—. (13) 



2 \ ^(^ + g-2) 2' 

From Eq.([13|) one can see that there are no solutions with arguments in the interval < 9 < 6 , 
where 8 Q = 2 arccos \J B~ x . This is the well known gap pi in the distribution of Yang-Lee zeros. 



It points to the absence of phase transitions in an ID ferromagnetic Potts model for Q > 1 at 
any real temperature. This is in good agreement with recent studies by Glumac and Uzelac [^4[ , 



and Kim and Creswick 0, where the Yang-Lee zeros of the ID Potts model was studied by 
using the transfer matrix method. Comparing our formula ( fL3|) with the corresponding formula 
(14) of Ref. H one can see that the argument of the derivative in our method is nothing but 
the difference in the arguments of two maximal eigenvalues in the transfer matrix method. It 
follows from the Eqs.flTTD and ( |T3"D that the Yang-Lee edge fields correspond to = and have 
the form 

1 



^± = - 2 \^(z-l)(z + Q-l)±^l-Q^ , (14) 

and /x± are complex for Q > 1. 

B < 1 when Q < 1 and it is positive or negative depending on z. Hence, B is negative when 
z < 2 — Q (Q < 1), and all values of /x's in ([TT|) are real and lie between the /x_ and /x + where 

fl± > 0. 

1 Zeros of the partition function considered as a function of complex Q. 
2 The results will be published elsewhere. 
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< B < 1 when z > 2 — Q (Q < 1), and all values of /x's in ( [LTD are either real or complex 
depending on 0. For < < 0o, where O = 2 arccosv^B, solutions ( [LTD are real and lie in 
the interval For O < < n the solutions (|TTJ) are complex conjugate and lie on the 

circle with radius R = z+ ®~ 2 , 



For Q > 1, the differentiation of both sides of Eq.(]T3T) with respect to /i and gives the 
density of Yang-Lee zeros g(8) in the form 

(6) ^ l s ^ n fl 

27r */sin 2 f -sin 2 ^ 



From Eq.(|l5D it follows that the density function g{6) for Q > 1 diverges in the gap points [i± 
with the exponent a = — |, i.e. g(#) oc \6 — 9 \ a when — > or — > 6* . 

For Q < 1, the corresponding density function may be obtained by differentiation of 
both sides of Eq.flnj) 

= — ; l^-v^'^ (16) 
2^ A /(, /+ - /i )(/ i -0_) 

One can see that g(/i) diverges in the points /i±, i.e. oc |/x — fi±\ a , with the exponent 

a = — |. Thus, for Q < 1 the density function of Yang-Lee zeros of the ID Q-state Potts 
ferromagnetic model has singularities only at the points fi± [Eq.([H])] with the edge singularity 
exponent a — — ~. The same is true for the antiferromagnetic case. 

Glumac and Uzelac |14|] considered for Q < 1 case the contribution of the third eigenvalue 
of the transfer matrix (A2 in their notations). We want to note that for the ID Potts model 
the transfer matrix method gives three and more eigenvalues only for Q > 2, hence, the 
third and other eigenvalues should be neglected for Q < 2. In this case the study of two 
maximal eigenvalues gives the same results as our method. The summary of results of the ID 
ferromagnetic Potts model is given in Fig. [| It is interesting to note that the argument <fi of 
the derivative in Eq. @ corresponds to the argument of the maximal eigenvalue in the transfer 
matrix method. 

The antiferromagnetic case may be studied in the same manner. The results are shown in 

Fig. a 



4 Yang-Lee zeros of the Q— state Potts model on a Bethe 
lattice (Q > 1) 

Let us, at first, consider the ferromagnetic Potts model on the Bethe lattice with coordina- 
tion number 7 = 3. In this case the system (H) may be studied analytically for neutral fixed 
points (k = 1). The exclusion of x from both equations (|5p gives the following equation 

z\Q - l)fx 2 - [2{Q - 1)2/(4 cos0 + 1) +y 2 (4e i * sin 2 f + 1) - 2{Q - l) 2 )fx + (z + Q - 2) 3 = 0. (17) 

where y = (z — l)(z + Q — l)/2 and <fi G [0,27r]. This equation describes the border of the 
metastability region in the complex \i plane (Fig. |J). The dashed areas in Fig. [| show the 
metastability regions for different temperatures. Inside the metastability region there are two 
attracting fixed points and there is only one outside of it. The other fixed points are repelling. 
Note that such a behavior is valid for any 7. At p,± at least two of the fixed points become 
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neutral. It will be shown below that at fi± the magnetization function ([?]) is singular and 
fi± correspond to the Yang-Lee edge singularity points. These are solutions to the equation 
(P~T|) for = and the critical temperature may be obtained from the condition /i + = It 



follows from Eq([T7|) that the edge singularity points lie on a circle with radius R 2 — ^ z+< ^ ^ 



V z 3 (Q-i) ' 

Since there is no phase transition on the boundary of a metastability region it will not give 
rise to zeros of the partition function P6| . The metastability region in the complex /z-plane 



points to the existence of the first order phase transition for complex magnetic fields. The 
problem of finding the Yang-Lee zeros of models with first order phase transitions attracted 
much attention for many years. Recently, Biskup et al., showed that the position of partition 



function zeros is related to the phase coexistence lines in the complex planes [231]. Dolan et 
al., used this approach to study Fisher zeros for Ising and Potts models on non-planar ('thin') 
regular random graphs. It is interesting to note that the locus of Fisher zeros on a Bethe lattice 
is identical to the corresponding random graph |27| . For our models the phase coexistence line 
is defined as a line in the complex plane, where the absolute values of the recursive function 
derivatives in two attracting fixed points are equal(see also |16[). Our numerical study showed 
that the phase coexistence line for T > T c is an arc of a circle with radius ending at the 
edge singularity points (the dashed line in Fig. |J). 

Let us now study the analytical properties of the magnetization function (|7|) to prove that 
the edge singularity points /i± correspond to its singularities. The fixed point equation of the 
Potts-Bethe mapping ([5]) may be written in the form 

^^.(q-p'-t'+q-g). (18) 

1 — Z X 

For continuity at /i ^ 0, x is defined to be equal to (Q — l)/(z + Q — 2) at /i = 0. Considering 
x as function of /i, one can study the singularities of x{fi) that also correspond to singularities 
of M(/x). The singular points of x(/i) are \i = 0, fi = oo and (i{x±), x± being the points, where 
the derivative dx/dfi is infinite. x± satisfies the equation 

^ + ^ + q-2)(2-7)-7(q-i) j + £ ^-2 = 

2(Q-1)(7-1) z(Q-l) 

Eq.fllDD is obtained after differentiation both sides of Eq.(fIB|) with respect to fi and x from 
the condition that dfi/dx vanishes. Note that Eq.([l9|) may be derived from the system (§) for 
k — 1 and = by excluding \x from both of the equations. At the singularity points of M(/x) 
the recursive function has two neutral fixed points with = and the others are repelling. 
The critical temperature may be obtained from the condition /i + = /i_ or setting the 



discriminant of the quadratic equation (19) to zero 



i(2-Q + ^(Q-2)2 + 4(Q-l) 7 2 /(7-2) 2 )) for Q > 1, 
l-Q for Q < 1. 

For Q > 1 and 1 < z < z c the solutions x± are complex conjugate numbers with modulus 
R x = \J z ^qZ^ and for z > z c they become real valued. For Q < 1 the solutions x± are always 
real numbers. Substituting x± = R x e lOL± into Eq.(|18|) one finds after some algebra 

H ± = Re i[S±+a±( -' r - 1) \ (21) 
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where 



R^RHQ-1), tan%= tan ^ and cosa ± = ^^T^ 2 - (22) 

Note that the fixed point equation f(x) = x and the magnetization M(fi) are invariant under 
the transformation G:{yU— >— , x —}. Moreover, for the magnetization function one has 




M(/i) = — M — • (23) 



The Yang-Lee zeros of the Potts model on the Bethe lattice may be obtained also analytically 
by studying the analytical properties of the magnetization function (0). This was done for the 



first time by Bessis, Drouffe and Moussa [24] for the Ising model on Bethe lattice. Making use 



of singularities of x(fi) in analogy with P4J a careful analysis of the analytic properties of the 



function M M Eq. ([?]) for Q > 1 gives the following picture for Yang-Lee zeros: (a) For z < z c : 
M(n) is analytic in the complex plane cut along an arc of the circle with radius R^ = Jy(g^_iy}-$ > 
that contains the point — R^ and is limited by the points n±, which are complex conjugate. 
Due to relation (|23|), the discontinuity along the cut is real and never vanishes except at the 
Yang-Lee edge singularity points This is in a good agreement with the Yang-Lee theory ||. 
The density of zeros g(Rfj,, 9) may be calculated from 

lim M(/i) ^ =re w - lim M(/i) = -4irg{R p ,6). (24) 



*Ru 



Our numerical study of Eq. (|24D shows that close to the edge singularity points fi± = e ±id °, 
gi^R^yO) has an exponential behavior g(R^,6) oc \9 — 9 \ a and Oq{T) ^ (T — T C ) A where a = 1/2 
for z < z c , a = 1/3 (cr = |) for z = z c and A = | (A = (35). Our results are in a good 
agreement with those of [J22| , [ |21| , where it was shown that the a exponent is universal and 
always equal to 1/2 for T > T c and a — 1/5 for T = T c in ferromagnetic models on lattices 
with spatial dimension d > 6. 

(b) For z > z c the function M(fi) is split into two different functions: M + (/j.) defined for 
|/i| < R^ and M_(/i) defined for \jj,\ > R^. The function M + (/i) can be analytically continued 
outside the circle = R^ into the z plane cut along the real axes from /i(x + ) to oo, where 
x + is the largest root of Eq.(|l^). The point fi{x + ) increases from unity to infinity when z 
increases from z c to infinity. The discontinuity of M + (/i) across the cut is purely imaginary 
and does not change the sign. Hence, it has no influence on physical properties of the model 
and the Yang-Lee zeros lie on the circle \^\ = R^. The points fi± correspond to the boundary 
of metastability for real magnetic fields and the plots of the iteration function f(x) at these 
points are given in Fig. |5|. 



5 Complex Magnetic Field metastability regions 

In the previous section the neutral fixed points of the recursive function f(x) was con- 
sidered only. It was found that the set of values of the magnetic field, for which the recursive 
function (^) has at least one neutral fixed point gives the boundary of the metastability region 
in the complex //-plane. Inside it there are two attracting fixed points and others are repelling. 
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Numerical experiments show that in the //-plane the recursive function Eq.([5]) has a complex 
behavior with period doubling bifurcations. The question arises: What will happen to the 
metastability region if neutral periodical points of period k > 1 are also considered? To answer 
this question one has to study the system (H) for any k > 1. Since it is impossible to solve 
the system (|8]) directly for k ^> 1 and 7 > 3 analytically and even numerically for large k and 
7, the method developed in |2(J is used. It gives a numerical algorithm for searching neutral 
periodical points for the recursive functions like Eq.(|5]) and is based on the theory of com- 
plex dynamical systems and the well known fact that the convergence of iterations to neutral 
periodical points is very weak and irregular, i.e. one has to make a number of iterations in 
order to approach a neutral periodic point. The algorithm is to find all critical points of the 
recursive function and investigate the convergence of all the orbits started at critical points 
(critical orbits). If all critical orbits converge to any attracting periodical point one says that 
the recursive function has only attracting and repelling periodical points. If at least one of the 
critical orbits does not converge, for example, after n iterations, one says that the recursive 
function has a neutral periodical point. Of course, the last statement is not rigorous from the 
strong mathematical point of view since a weak convergence to an attracting periodical point 
is also possible. Depending on the choice of n and e (the accuracy of approaching an attracting 
periodical point), the resulting picture on the //-plane may change. Our simulation experiments 
show that n = 700 and e = 10 -5 are optimal values and the data generated with this algorithm 
do not qualitatively change when n and/or e differ from their optimal values. For more details 



of the method and the C + + program code see [20, Eq 



One can easily find all critical points of the mapping / from Eq.(|^). The critical points are: 
x = with multiplicity 7 — 2 and x = 00 with multiplicity 7 — 2. The degree of our mapping 
/ is d — (7 — 1). It may be shown that these are the only critical points of the mapping 
/ [PU| . Hence, one has to consider only the orbits of the points Xq = and Xo = 00. For 
numerical calculations it is convenient to start iterations at the points x\ = f(0) = 1/z and 
Xl = f(oo) = (z + Q-2)/(Q-l). 

In Fig. the dynamics of metastability regions of the Q-state Potts model on the Bethe 
lattice with coordination number 7 = 3 and z = 3 is shown depending on Q. We have 
experimental evidence that in white regions all critical orbits converge. Figs. ^|(a)-(d) show the 
metastability regions for the case Q > 1 . It is seen that the sets of black points are similar to 
the boundary of the Mandelbrot set of the quadratic mapping z — > z 2 + c (Figs. |(c),(d)). This 
fact is known as the universality of the Mandelbrot set [p8 [. The metastability region of the 
Ising model (Q = 2) at the critical temperature is shown in Fig. [](c). It intersects the positive 
semi- axis at /1 = 1, which is an evidence of the existence of real temperature second order phase 
transition in conformity with exact calculations [I8|. The Q < 1 case is shown in Figs, ^(e)-(i) 



One can see that Fig. ^|(e) resembles a mirror reflection of the Mandelbrot set boundary of (cf. 
Fig. |6|(d)-(e)). It is noteworthy that for Q = 1 and 7 = 3 the Bethe-Potts mapping becomes a 
quadratic one |29j, and our numerical method fails ||20|| . By lowering Q the metastability regions 
become more and more complicated Figs. |](f)-(i). Note that the dynamics of the metastability 
regions remains the same if one fixes Q and changes the temperature. 
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6 Conclusions 



We showed that the stability analysis of attracting fixed points of the recurrence relation, i.e. 
the condition that the absolute values of derivatives of the governing recurrence relation in 
two attracting fixed points are equal at a phase transition, may be successfully applied to a 
study of zeros of the partition function. Note that in the one-dimensional case the condition 
of a phase transition is equivalent to the condition of existence of neutral fixed points. It will 
be interesting to check the Monroe's conjecture that at a phase transition the box counting 
dimension of the Julia set of the governing recurrence relation is a minimum for zeros of the 



partition function [^3| . We suppose that the box counting dimension of the Julia set of the 
recurrence relation should be a minimum for zeros of the partition function also. This may serve 
as a new criterion for studying zeros of the partition function for models on recursive lattices. 

In conclusion we observe that numerical methods proposed in this paper are generic and 
may be used for investigation of zeros of the partition function and metastability regions of 
other models on recursive lattices. 
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Figure 2: A schematic representation of the Yang-Lee zeros of ID ferromagnetic Potts model. 
Here R = Z+ Q~ 2 and /x± are defined in (Pj). 
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Figure 3: A schematic representation of the Yang-Lee zeros of ID antiferromagnetic Potts 
model. Here R = 2 ~ Q Z ~ Z and /x± are defined in (PI). 
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Figure 4: Metastability regions and Yang-Lee zeros of the Q-state Potts model on the Bethe 
lattice with coordination number 7 = 3 and Q = 2. The solid lines correspond to the boundary 
of metastability regions (gray filled areas). Dashed lines present an arc or circles of radius 
R 2 ^ = }*(^_~^]-2 and correspond to Yang-Lee zeros, (a) T > T c (z = 1.8), /j± are Yang-Lee edge 
singularity points; (b) T = T c critical point (z = z c = 3); (c) T < T c (z — 6), /jl± are spinodal 
points of the model. For more details see the text. 
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Figure 5: Plots of f(x) function ([5|) indicating the existence of neutral fixed points for different 
temperatures and magnetic fields at Q = 2 and 7 = 3. (a) T < T c (z = 6) and exp(h/kT) = 
(b) T = T c (z = 3) h = 0; (c) T < T c {z = 6) and exp(hjkT) = //+. 
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Figure 6: The dynamics of metastability regions of the Q-state Potts model on the Bethe lattice 
with coordination number 7 = 3 and z = 3 for different values of Q. For more details see the 
text. 
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